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ABSTRACT 

Wc present a new systematic analysis of the early radiation era solution in an in- 
teracting dark energy model to find the adiabatic initial conditions for the Boltzmann 
integration. In a model where the interaction is proportional to the dark matter den- 
sity, adiabatic initial conditions and viable cosmologies are possible if the early-time 
dark energy equation of state parameter is w e > —4/5. We find that when adiabatic- 
ity between cold dark matter, baryons, neutrinos and photons is demanded, the dark 
energy component satisfies automatically the adiabaticity condition. As supcrnovae la 
or baryon acoustic oscillation data require the recent-time equation of state parameter 
to be more negative, we consider a time- varying equation of state in our model. In a 
companion paper arXiv:0907.4987J wc apply the initial conditions derived here, and 
perform a full Monte Carlo Markov Chain likelihood analysis of this model. 

Key words: cosmology: theory, cosmic microwave background, cosmological param- 
eters, dark matter 



1 INTERACTING DARK ENERGY 

Dark energy and dark matter, the dominant sources in the 'standard' model for the evolution of the universe, are currently 
only detected via their gravitational effects. This implies an inevitable degeneracy between them. A dark sector interaction 
could thus be consistent with current observational constraints. We look at such a model, assuming that the dark matter and 
dark energy can be treated as fluids whose interaction is proportional to the dark matter density. 
For interacting dark energy, the energy balance equations in the background are 



Pe = 
Pde = 

where 7i 



-3Hp c + aQ c , (1) 
-3H(1 + w de )p de + aQ de , Qde = -Qc (2) 

a' /a, Wde — Pde/ Pde is the dark energy equation of state parameter, a prime indicates derivative with respect to 



Various forms for Q c have been investigated (see, e.g. Wetterich 


1995): Amendola (1 19991): Billvard & Colev (200C): 


Zimdahl & Pavonl d2001); Farrar & Peebles! d2004); Chimento et al. (200; 


): Olivares et al. d2005): Koivisto (20051): Guo et al. 


(2007 
(200fi 


); 
); 

): 


Sadiadi & Alimohammadi (hood): Boehmer et al. (2008): He & Wans (hoosl): Ouartin et al. ( 2008TI: IPereira & Jesus 


Ouercellini et al. (boost): IValiviita et all (hoOct): llle et al. d200Sh: 


Bean et al.l (20081): Chonechitnanl ()2009l): Corasaniti 


(200c 


Caldera-Cabral et al. J2OO9I): Gavela et al.l (20od): Jackson et al. 


(2009|)). We consider models where the interaction 



has the form of a decay of one species into another - as in simple models of re heating and of c u rvato n decay jMalik et al 



2003 



Sasaki et allbOQfj ; lAssadullahi et al.ll2007l ). Such a model was introduce d bvlBoehmer et al.l (|2008T ); I Valiviita et al.l (|2008T ). It 
is not derived from a Lagrangian [in contrast with e.g. Wetterichl ( 1995 ): Amendola ( 19991 )]. but it is motivated physically as 
a simple phenomenological model for decay of dark matter particles into dark energy. In this sense, it improves on most other 



phenomenological models, which are typically designed f or mathematical simp licity, rather than as models of interaction. The 
methods that we use here and in the companion paper ( Valiviita et al .1120091 ) may readily be extended to other interactions , 



including those based on a Lagrangian. We assume that in the background the interation takes the form (|Boehmer et al.ll20Q; 
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Valiviita et ailboosl ) 



Qc = -r Pc , (3) 

where T is the constant rate of transfer of dark matter density. Positive P corresponds to the decay of dark matter into dark 
energy, while negative r ind icates a transfer of energy from dark energy to dark matter. 

In lValiviita et al.l (|2008l ) we considered the case of fluid dark energy with a constant equation of state parameter — 1 < 
Wde ^ —4/5, and found a serious large-scale non-adiabatic instability in the early radiation era. This instability grows stronger 
as Wde approaches —1. Phantom models, Wde < — 1, do not suffer from this instability, but we consider them to be unphysical. 

The instability is determined by the early-time value of Wde ■ We will show that the models are not affected by the large- 
scale non-adiabatic instability during early radiation domination if at early times Wde > —4/5. If we allow Wde to vary, such 
a large early Wd e can be still consistent with Supernovae Type la (SN) and baryon acoustic oscillation (BAO) observations, 
provided that at late times Wde ~ — 1- In this paper, we represent Wde via the parametrization Wde = wo + w a (1 — a) 
( Chevallier fc PolarskilbOQll : lLinderll2003h . which we rewrite as 

Wde — w a + w e (l — a) , (4) 

where w e = wo + w a is the early-time value of Wde , while wo is the late-time value. 

There are two critical features of the analysis of interacting models, which are not always properly accounted for in the 
literature: 

• The background energy transfer rate Q c does not in itself determine the interaction in the perturbed universe. One should 
also specify the momentum transfer rate, preferably via a physical assumption. We make the physical assumption that the 
momentum tra nsfer vanishes in the da r k matter rest-frame; this requires that the energy-momentum transfer rate is given 



covariantly by ([Kodama fc Sa saki 1984; IValiviita et al.ll200 



(5) 



Q c = -r Pc (i + »5 c ), 

where u% is the dark matter 4- velocity, and <5 C = Sp c /p c is the cold dark matter (CDM) density contrast. 

• Adiabatic initial conditions in the presence of a dark sector interaction require a very careful analysis of the early-radiation 
solution, both in th e background and in the perturbations. We derive these initial conditions by genera lizing the methods of 
Doran et all |2003h to the interacting case, thereby extending our previous results (|Valiviita et al.ll2008h . 



We give here the first systematic analysis of the initial conditions for perturbations in the interacting model given by 
Eq. ( [ 5)l - and our methods can be adjusted to deal with other forms of interaction. In the companion paper IValiviita et al. 
( 2009;) we report the results of our full Monte Carlo Markov Chain likelihood scans fo r this mod e l. Cosmological perturbation s 



of other interacting models h ave been investigated, e.g., in 



Jackson et al.1 j200gF: 
( 2009h ; lKovama et al 



Gavela et al. "(120091); lLa Vacca et al 
(|2009h ; iKristiansen et al.l |2009l) 



Amendola et al 



Mainini fc Bonomettol ]2007l);lBean et al.1 (120081); I Vergani et al 



ihOQSh: ISchaferi tOQch: ISchaefer et al] (120081): iHe et al.1 faood): Bean et al.1 (hoosl): ICorasanitil ihoosh: Ichongchitnanl ]2009l) 



oogfk IE 



( 20031); iKoivistol J2OO5I) ; lOlivares et~al] tod) 



39) 



2009l);|Pettorino fc Baccigalupil j20oi);lLa Vacca fc Colombo! 



(2009); IHe et al 



2 PERTURBATION EQUATIONS 

The scalar perturbations of the spatially flat Friedmann-Robertson- Walker metric are given by 

ds 2 =a 2 | - (l + 2(t>)dT 2 + 2d l BdTdx l + [(1 - 2iP)8 tJ •+ 2d i d j E^ dx'dx 3 1 . (6) 

In the perturbed universe, the dark sector interaction involves a transfer of momentum as well as energy. The covariant form 
of energy-momentum transfer for a fluid component A is X7„T^ V = Q^, where = a~ 1 (Q e ,0) = — Q% in the background. 
The perturbed energy- momentum transfer 4- vector can be split as ( Valiviita et al. 20081 ) 



tf = -a[Q A (l + <t>)+5Q A ], 



h = ad t ( f A - Qa-^ 



(7) 



where k is the comoving wavenumber, /a is the intrinsic momentum transfer potential and 9 = (p + p) _1 XXpa + Pa)Qa is 
the total velocity p erturbation (6 = — k 2 v). The evolution equat ions for density perturbations and velocity perturbations for 
a generic fluid are ( Valiviita et al. 2008; K odama fc Sasaki 1984 ) 



8' A + 3H(c 2 A - w A )S A + (1 + wa)0a + 3H[3H(1 + w A )(c 2 A - w A ) + w' A ] 

clQa 



3(1 + w A )i>' + (1 + w A )k 2 (B - E') 



9' A +H(\-ic 2 sA )e A 



(1 + w A ) 



k 2 S A + 



PA 

2w A 



S A + 3H(c 2 A - w A 



3(1 + wa 



■ k 2 Tv A — k 2 (f> 



k' 2 
aQ A 



a 

H 

PA 



(1 + W A )p A 



1 + cJa) 



(1 + w A )p A 



(8) 



k 2 f A .(9) 
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Here c 2 A is the sound speed, and tta is the anisotropic stress. For our model -Kde = 0, and w e set c 2 de = 1, as in standard 
non- interacting quintessence models, in order to avoid adiabatic instabilities (see discussion in Valiviita et al. ( 20081 )). 
For the interaction defined by Eq. ((5j), we find from Eq. © that 



f c = rg (0c - 6) = -f de . 

Then we can write the dark energy and cold dark matter perturbation equations for our model: 

S' de + 3H(1 - w de )8 de + (1 + w de ) [9 de + k 2 (B - E')] + 9H 2 {1 - wle)^ - 3aH 2 w a ^- ~ 3(1 + w de )^' = 

Pc 



aF- 

Pde 

— 2H.6 dl 



S de + 3H(1 - Wde) + ' 



k 



(1 + w de ) 
5' c + e c + k 2 (B-E')-3ip' 
6' C + H0 C - k 2 4> = 0. 



8 de — k <f> = 



aT 



pc 



(1 + W de ) pde 



— 28de 



-aTtf) , 



(10) 



(11) 

(12) 

(13) 
(14) 



3 BACKGROUND SOLUTION IN EARLY RADIATION ERA 



The background solution in the early radiation era (ptot — pr) is important for finding the initial conditions for the integration 
of cosmological perturbations. In what follows we use occasionally the Hubble parameter H = a" 1 ?^ instead of the conformal 
Hubble parameter TL. In the radiation era we have 

Tt = t -1 and a = HoV^rt) t , (15) 

pro/Pcvito ~ 2.47 x 10 _5 /i -2 is the radiation energy density 



where Ho is the conformal Hubble parameter today, and fi, o 
parameter today. Here h is defined by Ho = h X lOOkrns - Mpc" 1 , and as do = 1, we have Ho 
H = {2t)~ x and a = (2_ffoV^rO t) 1//2 where t is the cosmic time. By Eq. (|15|l we find 

t = {HoVfU /2)t 2 . 



Ho- Furthermore, we have 
(16) 



We define the ratio of dark energy to cold dark matter density r — pd e / p c - Then, employing Eqs. |T} and ([2]), 

1/21 

■ 

r + iw a 



3w e 
~2t 



Ho V ^rO 



2t 



r + T, 



(17) 



where the dot indicates derivative with respect to cosmic time. At early enough times, |r + 3u> a (-Ho\/^ro/2t) 1/ ' 2 | <S 3\w e \/(2t), 
and we can neglect the term in square brackets, so that the solution is 



t 



-3u) e /2 



2F 



+ 3^T2^ (18) 
where r ro f is an integration constant corresponding to p d e/pc at some (early) reference time t Te f in the case where F = 0. 
From Eq. (|18[) we find that we have two regimes, depending on the value of the early-time equatio n of state parameter w e . If 

We ^ 



2/3, then the second term dominates over the first as t — > 0, and we recover the solution of Valiviita et al. ( 200Sh : 



Pde 
Pe 



aT 



:H~ 



3w e + 2 

If We > —2/3, then the first term in Eq. (|18[) dominates: 

-3w a /2 



C(kr) 



Pde 
Pe 



r rB i 



t 
tret 



C{kr) 



G 



C 



HqT vn7o 
k 2 3w e +2 

Ho^/^lrO 

1 2t vc( k 2 



-3iu e /2 



For the background evolution of p c in the radiation dominated era, 



Pc 



27 + r )^ 



(19) 



(20) 



(21) 



the second term in brackets is negligible relative to the first at times t 3/(2F), or t < r sw i tc h = (\/^ro'|r / Hol/S) -1 ^ 2 Hq 1 . 
For these times, p c oc a -3 . In typical models that provide a good fit to Cosmic Microwave Background (CMB) data, H^ 1 = 
O(10) Gpc, and the conformal time at matter-radiation equality is r eq = 0(100) Mpc. If we demand that the evolution of p c 
is effectively standard during the whole radiation dominated era, i.e. r sw i tc h > r eq , we require 



Ho 



30 000 

V^rO 



W e h, 



(22) 



where h ~ 0.7. As we study in this paper coupling strengths |r///o| ^ 1, we can safely assume that the cold dark matter 
evolution during radiation domination is completely the standard non-interacting one p c = p% q (a/a eq )~ 3 , where p% q and a eq 
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are the dark matter energy density and the scale factor at matter-radiation equality, respectively. Noticing that the radiation 
energy density can be written as p r = p r q (a/a eq )~ 4 , and that p eq = p Eq by definition, we find that in the radiation dominated 
era 

p c a Ho vfirt) /ooN 
— = =LJ 2 kT, u 2 = -r , (23) 

pr Q*eq <^ CLeq 

where we used Eq. (|15[) . In the non-interacting case we could continue by setting a eq = i} r o/il c o, but in the interacting case 
the dark matter evolution from r sw i tc h up to today (to) differs from oc a~ 3 : by Eq. (|21|l . it follows that at recent times, for a 
positive r, the dark matter density decreases faster, and with a negative T it decreases slower than a~ 3 . Therefore we cannot 
do the "a -3 scaling" all the way up to today, but instead have to stop at some early enough reference time. Here we choose 
the time of matter-radiation equality. 

An upper limit to the early dark energy (DE) equation of state w e could be set by requiring dark matter domination 
over DE at early times. Then Eq. (|20|l would set the constraint w e < 0. However, if the DE equation of state is close to 
at early times, it could well mimic the behaviour of cold dark matter. On the other hand, if w e is close to 1/3, the "DE" 
component would behave like radiation at early times. So, for ^ w e < 1/3, we conclude that the fluid which at late times 
behaves like dark energy, behaves at early times like a combination of matter and radiation. As this case cannot be ruled out, 
we set a conservative upper bound on w e by demanding that in the early universe DE does not dominate over radiation, i.e., 
for t — > 0, we have pde/pr — * 0. Using Eqs. ([20]) and ([23)) , we find 

Pde = Pde_ Pc =C r wa ( feT )i-s«. (24) 

Pr Pc Pr 

which implies, as expected, w e < 1/3. 



4 SUPER-HORIZON INITIAL CONDITIONS FOR PERTURBATIONS 



In order to solve numerically the perturbation equations we need to specify initial conditions in the early radiation era. 
The wavelength of the relevant fluctuations is far outside the horizon during this period: fcr <C 1. To compute the initial 
conditions we start by writing the perturb ation equations of each species and the perturbed Einstein equations in terms of 
the gauge-invariant variables developed by Bardeenl (|l980l ): 



$ = -iP + H(B- E'), 
A A = S A +H- 1 ^, 

PA 



* = <t> + H(B-E') + (B- E')' , 
V A = fc _1 0A + k(B - E') , U A 



TV A- 



(25) 



The general evolution equations for the density, velocity and anisotropic stre s s pert urbations A a, V a and Yl A and the Einstein 
equations for the metric perturbations <E>, ^ are given in Kodama fc Sasakil ( 1984T ). 

We use and generalize the systematic method of Doran et al.l ( 2003 ) in order to analyze the initial conditions in the 
interacting DE model with time- varying w e . The results are derived below, but let us summarize the key points before going 
to the details. The conclusion is that we can use adiabatic initial conditions for 

(26) 



4 , 

— ^ We 
5 



2 

— < W e 

3 



1 

<3 



For both of these intervals, the initial conditions for all non-dark energy quantities are the same as in the non-interacting case. 
For the second w e interval, the initial dark energy density perturbation is the same as the standard one, Ade = 3(1 + w e )A 7 /4, 
whereas for the first interval, we find a non-standard initial condition Ade = A 7 /4. The difference arises because of the different 
background evolution in the two cases (as given in the previous section). Note that for w e < — 1 it is also possible to have 
adiabatic initial conditions, bu t we c onsider this case to be unphysical. For —1 ^ w e ^ —4/5, we recover the non-adiabatic 
blow-up case of lValiviita et al. |200Sl) . 

Similar considerations could be extended to the early matter dominated era. The key difference there is that the back- 
ground behaves differently for the interval —1/2 < w e < 1/3 than for w e ^ —1/2, where the interaction modifies the DE 
evolution. In the matter era, a non-adiabatic blow-up may thus happen if 



1< w e < -1/2. 



(27) 



A detailed analysis shows however that in the interval 



2/3 < We < -1/2, 



(28) 

the "blow-up" mode is in fact a decaying mode, and hence (non-standard) adiabatic evolution on super-Hubble scales is 
possible. In the interval —1 < w e < —2/3 the non-adiabatic "blow-up" mode is rapidly increasing, and will dominate unless 
|T| is suitably small. Therefore, a blow-up in the matter era will make large interaction models with —1 < w e < —2/3 non- 
viable, while the blow-up in the radiation era ruins all interacting models (no matter how weak) with —1 < w e < —4/5. We 
summarize these results in Table [1] 
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Table 1 . The evolution of perturbations on super-Hubble scales with various values of the early dark energy equation of state parameter 
in the radiation and matter dominated eras (RD and MD respectively). "Adiabatic" means that it is possible to specify adiabatic initial 
conditions so that the total gauge invariant curvature perturbation £ stays constant on super-Hubble scales. "Adiabatic (standard)" 
means that the behaviour of perturbations at early times on super-Hubble scales is the same as in the non-interacting model. 



Wde. in the RD or MD era Radiation dominated era (RD) Matter dominated era (MD) Viable? 



w de < -1 
-1 < w de < -4/5 
-4/5 Sj w de < -2/3 
-2/3 Sj w de < -1/2 
-1/2 < w de < +1/3 



adiabatic 

"blow-up" isocurvature growth 
adiabatic 

adiabatic (standard) 
adiabatic (standard) 



adiabatic 

"blow-up" isocurvature growth 
isocurvature growth 
adiabatic 

adiabatic (standard) 



viable, but phantom 
non-viable 

viable, if T small enough 

viable 

viable 



Assuming tight coupling between photons and baryons, so that Vt = Vy, passing from conformal time r to the tim e 
variable x = kr, using a rescaled velocity Va = Va/x and a rescaled anisotropic stress ft a = Ha/x 2 , as in iDoran et al"1(|2003r i. 
we obtain the following evolution equations: 

dA c 



din a; 

dV c 
din a; 
dA 7 
din a; 
dVy 
din a; 
dA b 
din a: 
dA„ 
din a; 
dV v 
din a; 
rfft„ 
din a; 
dA de 
din a; 



dVde 

din a; 
Here a 

$ 

d$ 
din a: 



"Ho 



= -xVo-— ax 2 \^-{l + w)V + 2U v Yl v + 2^ 



= -2V C + * , 
4 2f/ 

= 7A 7 -K,|flJ„ + 2f ; 



-x 2 V, 



4 2 f/ 

-A v -%- \x 2 tl v + Q„ft„ + 2* . 
4 6 



\v v -m v , 

5 

3(u> e - 1) { A de + 3(1 + w e ) [-* - SlvfiJ + (1 



3(We - 1) 



V,, 



, r 2 pc 

+ TT ax — 

rlO Pde 



A c — A x + 3(1 — WejVde ~ {3w e 



5) [* + fi„ft„] +^(l + w)V 



+ V de + 3Q„TL V +4* + -Laa, 2 -^ 

L + W e Tio p de 



2V de - n v n v 



1 + W e 



(Tio/k) 2 \/Q r (,, and we used the Einstein equations 

3 -2 

2" 

3 -( 
2 V 

-* - 3wU. 



x~ 2 |A + 3(l + w) [V"-*]} , 
l + w)V, 



to eliminate <3? in favour of *]/. 

Using Eqs. (J39J) and (gTJl, we have 

* = - 



E^, t , 7 ,», fe "a [Aa + 3(1+ wa)Vs 

J2A= c , b ,^, d£ Hi + ^A)n A + lx 2 



- fi„iu 



The total velocity appearing in Eqs. (|29[1 and (|37(l is 

V(l + u;)= Y, n A (l + w A )V A . 

A=c,b, r y,i/,de 



(29) 

(30) 
(31) 
(32) 
(33) 
(34) 
(35) 
(36) 

(37) 
(38) 

(39) 

(40) 
(41) 

(42) 
(43) 



Recalling that p = Pa, (1+w) = J2^a(1 + w a ), A = ^OaAa, and ft = J] ^aAa = fi„ft„, we then see that Eqs. (f29ti38]l 
form a set of 10 linear differential equations for 10 perturbation variables Aa, Va and ft„. (Note that IIa = for A 7^ v.) 
Since we are interested in the early radiation era, we make the approximation Q,a = Pa/p — PA/pr- Using Eqs. (|15[) , 
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(|19[) . (|23[) and l)24[l. and the (standard non- interacting) background evolution of photons, baryons and neutrinos, we obtain 

fi c — UJ2 x , 



fide = CtJ2X for W e ^ —2/3, 



for -2/3 < w e < 1/3, 



and fide = Cui2X 

fii/ = P„/Pr = , fi 7 = 1 - fib - fic - fide - fii/ ■ (44) 

The next step of the method proposed in iDoran et al. (2003) consists in writing the system of differential equations 



2911381) in a matrix form: 



dU 
din a; 



= A(x)U(x) , 



(45) 



where 

U T = {A c , Vc, A 7 , %, A b , A„, %, fl,, A de , V de j , (46) 

and the matrix A (a;) can be read from Eqs. (|29H38|I after substituting Eqs. (|42H44[) and the background evolution of p c /pde 
from (I19|) or (|20|) . depending on the value of w e . 

The initial conditions are specified for modes well outside the horizon, i.e. for i« 1. There will be several independent 
solution vectors to Eg. 1451 that we write as x x 'U^K If no term of A(x) diverges for x — > 0, then we can approximate A by a 
constant matrix Ao = lim^^o A(x). If we require more accuracy we can expand A(k) up to a desired order in x. For example, 
up to order x 3 the matrix A (a;) contains the constant term Ao as well as terms proportional to x, x 2 , and x 3 in the case where 
We ^ —2/3. However, in the case —2/3 < w e 1/3, A(a;) contains in addition to integer powers of x also non-integer powers 
1 — 3w e , 2 — 3w e , 3 — 3w e , etc., and 2 + 3w e , 3 + 3w e , 4 + 3w e , etc. The listed ones and possibly their multiples can fall in the 
range (0,3). For a given w e , however, one should drop those which turn out to be higher order than x 3 . 

Thus going beyond zeroth order, up to order x 3 , we can expand A and each solution x Xi U^ as 



A(x) 

c x 'U (l) (x) 



A + Ai x + A 2 x 2 + A 3 x 3 + 



j = ln-1 

E7$° + U[ l) x + U ( 2 l) x 2 + C/.y 



3=0 



r(0 „™(l-3^e)+j 

Bnj 1 



(i) 2+3m e +j 



(47) 
(48) 



AoC/q 



c/r = - 



U 



(i) 



1/ 



(i) 



where Aj, B nj , and Cj are constant (not depending on the time variable a;) matrices, and Uj are constant vectors. Note 
that for w e ^ —2/3 all B ?lJ , Cj, U^ n j, and Uq^ terms vanish. For simplicity, we demonstrate below this case, which leads to 
only integer powers. Substituting Eqs. (|47[) and (|48[) into the evolution equation (|45|l and equating order by order, we obtain 

'Jq , i.e. A; is an eigenvalue of Ao, and Uq is an eigenvector of Ao, (49) 
Ao - (Ai + I)!]" 1 [AiC/W + AoC/f] , (50) 

Ao - (Ai + 2)!]- 1 [a 2 C7 ( ( ; ) + ArC/f + A £/ 2 l) ] , (51) 

Ao - (Ai + 3)!]- 1 [a 3 C/W + A 2 [/$ l) + Ail/« + A C/^] . (52) 

Now the general solution to the differential equation (|45[) is a linear combination of solutions x x 'U^: 

U(x)=^2ax Xi U (i) (x), (53) 

i 

where Cj are dimensionless constants. If we define an initial reference time fmit, then the constants Cj = CiX X ^ it represent 
the initial contribution of the vector U^' to the total perturbation vector U{xi n it). The imaginary part of Ai represents 
oscillations of x Xi U (x), while the real part gives its power-law behaviour: x Xi UQ\x) = a; Ro(A * ) cos[Im(A i ) lna;]. The 
contribution corresponding to the eigenvalue(s) with largest real part, Re(Ai), will dominate as time goes by, while initial 
contributions from eigenvectors corresponding to Ai with smaller real part will become negligible compared to the dominant 
mode. Hence, to set initial conditions deep in the radiation era but well after inflation, it is sufficient to specify the contribution 
coming from mode(s) with largest Re(Ai). 

From now on we divide the treatment into two cases — 2/3 < w e < 1/3 and w e ^ —2/3. Before proceeding, we should 



point out that the matrix method presented in lDoran et al.l (|2003l ) and applied to non- interacting constant-wide dark energy, 
represents a systematic and efficient approach for finding initial conditions. Once the matrix A(k) has been read from the set of 
first order differential equations (|29H38[) , one can feed it into a symbolic mathematical program such as Maple or Mathematica 
and easily extract the constant part Ao as well as the other parts (such as Ai, A 2 , B n j, Cj) up to any desired order. Then 
it is simple linear algebra to find the eigenvalues Ai and eigenvectors Uq of Ao and, if higher order solutions in fcr are needed, 
to substitute these step by step into Eqs. (|50[) - (|52p etc. in order to find the solutions x Xi U^\x). 
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4.1 Case -2/3 < w e < 1/3 



We substitute from Eq. (|42p . V from Eq. (|43jl . and the energy density parameters from Eq. (|44p into Eqs. 
using Eq. (|20[) for p c / pde in the last two of them, and taking the limit x — > 0, we find the Aq matrix: 



Then, 



/ 

-2 















\ 

where N = R v — 1. 

The eigenvalues of Aq are 
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(54) 



Ai = 



-2,-1,0,0,0,0,-1 -^HfM.-l + ^Ef^^-^ 



where 



-2 + 3to g . V-20 + 12iu e + 1 



(55) 



(56) 



2 2 

For the range < Rv < 0.405 and —2/3 < w e < 1/3, all eigenvalues have a non-positive real part. In (|56[) the term inside the 
square root, —20 + 12 w e + 9 w e 2 , falls between —24 and —15, and hence Re(A^) = — 1 + 3 w e /2, falls between —2 and —1/2. 

As explained in iDoran et all (120031 ). since the eigenvalue with largest real part, A = 0, is fourfold degenerate, it is possible 
to choose a basis from the subspace spanned by the eigenvectors with eigenvalue A = 0, so that physically meaningful choices 
can be made for the initial condition vector. One can form 4 independent linear combinations from the four vectors with A = 0. 
The physical choices are an adiabatic mode and 3 isocurvature modes. Here we choose adiabatic initial conditio ns, specified 



by th e condition that the gauge-invariant entropy perturbations Sab of every A, B =7, v, c, b vanish, where (jMalik et al 
2003! ) 



Sab = -3H^A A + 3H^-A B . 
Pa Pb 



(57) 



We will show later the interesting new result that demanding adiabaticity between the standard constituents automatically 
guarantees adiabaticity with respect to DE. 

We should remind the reader that for the interacting constituents the coupling appears in the continuity equation, and 
we should not use blindly the standard result 



Sab ~ 



(58) 



(l + w A ) (1 + w b )' 

where the 1 + wa factors result from applying the continuity equation to p' A / pA- Indeed, for cold dark matter in the early 
radiation era we find, using Eqs.JT}, ((3]), and (|15|l . 



- 3H^A C = ^ , 

P'c (1 + Wo) + (tf \/fU/3fc 2 )(fcr) 2 



where w c = 0. For DE, we find using Eqs.©, J3j, (fl5]l. and (gOJ 



3H^A de = 

Pde 



A de 



(1 + Wde) - (HoV^U)/3k 2 C)(kry+^ 



(59) 



(60) 
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At zeroth order in x = fcr we incidentally regain the standard non-interacting result (|58l) . From Sab 



it then follows that 
(61) 

Imposing this condition on a linear combination of the four eigenvectors with eigenvalue A = we obtain 

\ 



3 3 

Ac = A b = -A 7 = -A„ . 

4 4 



rr(adi) 
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-(5/4)7? 
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= Ci 
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A de 










(3/4) (1 + w e 
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Vdc 




adiabatic 






-(5/4)7? 



(62) 



/ 

where V = (15 + AR V )~ , and Ci is a dimensionless normalization constant corresponding to, e.g., A 7 and A„ at the initial 
time. 

I ■ vector (|62|l is identical to the standard adiabatic initial condition vector (see Doran et al. ( 20031 )1. In particular, it 
should be noticed that although we did no t require adiabaticity of DE, (|62p automatically satisfies the condition Sde,A — 
for all A = 7, v, c, b. In Doran et al. (2003)) this result was found for non-interacting dark energy. Here we have now shown 
that also interacting dark energy is automatically adiabatic, once cold dark matter, baryons, photons and neutrinos are set 
to be adiabatic. 

Finally, since all components of the vector (|62[) are different from zero, it is not necessary to compute terms to higher 
order in x, and we use Eq. (|62|l as our adiabatic initial condition for the computation of the CMB power spectrum for models 
wit h -2/3 < w e < 1 /3. 

Lee et al. ( 2006h have reported that the quintessence isocurvature mode decays away (in an interacting quintessence model 



which is quite similar to our set-up). After our systematic derivation of initial conditions, this decay can be tracked down to the 
fact that Re(Aj ) in Eq. (|55[) are negative. The reason for this i s that in quintess ence models the early-time equation of state 
parameter is typically larger than —2/3, indeed positive [but in Lee et all 1 20061 ) less than +2/3]. So Re(A^ ) is negative, and 
hence the isocurvature mode decays. In the next subsection we will see that also in the range —4/5 ^ w e ^ —2/3 (or w e < — 1) 
the DE isocurvature mode decays (although in this case the interaction affects the evolution of DE perturba tions), whereas 
in the range —1 < w e < —4/5 the DE isocurvature mode is a rapidly growing mode as recently realised by IValiviita et al 
( 20081 ). 



4.2 Case w e -2/3 

Si nce at early tim e s the equation of state can be approximated as a constant, w e = wo + w a , this case has already been studied 



Valiviita et alj ( 20081 ). where a constant —1 < Wd e ^ —2/3 was analysed. A serious non-adiabatic large-scale instability 



that excludes these models was found whenever —1 < Wde < —4/5, no matter how weak the interaction was. However, we 
notice that there is a limited region of parameter space, —4/5 ^ Wde ^ —2/3, where the instability can possibly be avoided. 
In the case of a constant DE equation of state parameter this range would be observationally disfavoured, since for example 
supernova data require that vide is closer to —1 at recent times. In the case of time varying Wde(a) we do not have this 
pr oblem as wp can be cl ose to —1 while —4/5 < w e ^ —2/ 3. In the following w e repeat the analysis of initial conditions done 



m 



Valiviita et al. (|200st ). but using the matrix method of iDoran et al.1 |2003h , extended to include the interaction, and give 



the conditions for a viable cosmology. 

Substituting >]/ from Eq. ((42}, V from Eq. (|43|l . and the energy density parameters from Eq. ((44} into Eqs. (f29H38|l . as 
well as pc/pde from Eq. (|19[1 into the last two of them, and taking the limit x — > 0, we find the Ao matrix, which is very 
similar to our previous result, Eq. (|54|l . This happens because everything remains unchanged, except that we must replace 
in Eqs. 1)370 and ((38} the evolution of p c /pde with Eq. (|19)l . and whenever Qde appears we must now substitute the oc x 3 
behaviour from ((44}, instead of the oc x 1 ~ 3w " behaviour. Therefore only the last two rows in (|54)l are modified, and will now 
read 

, + 3We q A/Xl+9^ 3A /-( u;e _ 1 ) - fl - (1 + 9 "' e) -3iMw e -l) -5 3(w e -l) 



c, 2+3m e Af(2+w s ) Af(2+w c ) c, fl„(2+m e ) R v (2+w c ) _ p 1 3+5m e 

U 4(1 + 1U C ) 1 + W e 4(l+1D e ) l + ™<! 1 + lUe l + W c 



(63) 



The eigenvalues of Ao are 
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where 



-5w e - 4 ± V3wl - 2 



± — ^V^e t (65) 

The first eight eigenvalues coincide with the previous case, Eq. (|55|l . Of those, four have a negative real part, corresponding 
thus to modes that will decay away quickly, and that we can neglect. The last two eigenvalues, Ajr , are instead very different 
from the previous case, and depend on the value of w e . The eigenvalue with a largest real part, Xg, is real and positive 
for — f < w e ^ — \/2/3. In addition to this, Re(A^) is positive also in the small r ange — \/2/3 < w e < —4/5. Therefore 
Re(A^) > for — 1 < w e < —4/5. This corresponds to the blow-up solution found in lValiviita et al. A+ is larger, the 

closer w e is to —1. There is no blow-up of perturbations for —4/5 ^ w e ^ —2/3, because then the largest Re(Ai) are zero. 



4-2.1 Case —1 < w e < —4/5; non-adiabatic blow-up 

We now calculate the initial condition vector U^ s '(x) corresponding to the fastest growing mode, Xg . At zeroth order in x, it 
is given by 

U { s)T (x) = {o, 0, 0, 0, 0, 0, 0, 0, -1 + s/Zwl - 2, l} . 



(66) 



In this case, since only the last two components of the vector are different from zero, we need to compute higher order 
corrections. It turns out that an expansion up to x 3 is necessary, and as explained both before and after Eqs. (|47[) and (I48|) , 
the expansion contains only integer powers of x, when w e ^ —2/3. Therefore we have 



A(x) 

U {9) {x) 



Ao + Ai x + A 2 x 2 + A 3 x 3 , 



V 



(s) 



U{ B) x + U\ g > x* + U 



(9) 



By substituting A(x) and x « U^ 9 '(x) into the evolution equation (|45|l and equating order by order, we obtain 



(;/) 



1)1] "AiC/o, 



(9) 



= - [Ao - (A+ + 2)1] 



U 



Ui s) = -[Ao-(A+ + 3)l] 1 



(A 2 U"o + AiE/i) , 
(A 3 Uo + Aat/i + AiC7 2 ) . 



(67) 
(68) 

(69) 
(70) 
(71) 



Using these formulas we find corrections to Eq. (|66|) . Keeping for each perturbation variable only the leading order (in x) 
terms, we obtain the following initial condition vector: 



C7 (g) (a 



Vc 

V, 

A,, 

A, 

V v 

11, 

A de 

V de 
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= Gi 





7 EjU3(l + u. c ) 3 [8fl„(l + mc )S-15(i 
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4(w e -l)(2+3w e )X 



-\-w 
2(2+3w e )M 




-15J-u 3 (l+u.e) 3 (2+2)i 3 
2(2+3uj c )A1 
4^0)3 (l+-ii. c ) 3 [3(7+8m, ! +-iiie 2 ) + (5+3ui e )S]x 3 

(2+ 3 UJ e )(S-2lDe)S 

Q 
1 



(72) 



/ 



whereas = w 2 (#oA) VHTo" = (H /k) 3 n r0 /a eq , Q, M and B are Q = V3w e 2 - 2-1, X = 5 [6 + 7w e + 3w e 3 + (3 + 5w e )Q\- 



4R V (1 + w e y (Q - 1 - 3w e ) andS = 8i?„ (1 + 



5w e + 3Q)]}.This 



_) 2 (5 + 3w e + 2g)+ 5 {9Q - 3 + w e [13 + UQ + 3w e (13 ■ 
solution coincides with equations (63)-(70) of Valiviita et al. ( 20081 ) . after substituting — A+ + 3, J = 1 — l6R v [5(n^ + 
2)(n^> + 1) + SR V ]~ , converting into Newtonian gauge (by using Eqs. (|25[) with B = E = 0) and conveniently renormalizing 
the vector. Equation (|72|l is the initial condition vector for the case —1 < w e ^ —4/5, when the dominant eigenvector is that 
corresponding to AjJ~. 

The initial condition (|72p is trivially adiabatic with respect to 7, v, c and 6, but not with respect to DE. Indeed, for DE 
we find using Eqs.©, ©, and (JT5J 

A de 
— 

Pde 

Therefore 



3H^A de 



(1+Wde) - (3w de + 2)/3 



3A fl 



(73) 



5',.; 



3A d 



C 2 Qx x i 



(74) 
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for any A = 7, v, c or b. Even if we were able to set the initial conditions at r = and demanded adiabaticit y there, after a 



short time the solution would not be adiabatic. Thus Eq. ([720 represents the non-adiabatic "blow-up" solution (|Valiviita et al 
2008) for the case -1 < w e < -4/5. 



4-2.2 Cases w e < —1 or —4/5 ^ w e ^ —2/3; adiabatic initial conditions 



In the range —4/5 < 



-2/3, as well as for w e < —1, we have Re(Ag ) < 0, so that the largest eigenvalue is the 



fourfold degenerate A = 0. [If w e — —4/5, then A = is fourfold degenerate, and there are also two oscillating solutions with 
Re(A^) = 0.] We look for a linear combination of the four eigenvectors (corresponding to A = 0) that satisfies adiabaticity 
(see Eq. (|6ip ~l of photons, neutrinos, baryons and cold dark matter. The resulting eigenvector is 



U 



(adi) 



Vc 

A- 
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A h 
A, 

Yr 
ft, 

A de 
V de 



\ 
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3/4 \ 
-(5/4)7> 
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-(5/4)7> 
3/4 
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-(5/4)7> 
—V 
1/4 
-(5/4)7> ) 



(75) 



adiabatic 

where V = (15 + 4ii„) _1 . This corresponds to equations (59-61) of Valiviita et al. ( 200S). All components e xcept 
equal to the initial conditions for —2/3 < w e < 1/3, Eq. (|62|l . However, as pointed out in lValiviita et al. I |200Sl ). A 

de 



corresponds exactly to the adiabaticity condition for DE: SdeA 
we find 



A de are 
= A 7 /4 

0. Namely, substituting the result (|73p into definition l|57|). 



Sde 



3A d 



3 , 
4^ 



(76) 



Thus Eq. (|75|) is an adiabatic initial condition vector for the cases —4/5 ^ w e ^ —2/3 or w e < — 1. 



5 CONCLUSION 



We have presented, for the first time, a systematic derivation of initial conditions for perturbations in a model of interacting 
dark matter - dark energy fluids, in the early radiation era. These initial conditions are essential for studying the further 
evolution of perturbations up to today's observables. They are the initial values for perturbations in any Boltzmann integrator 
which solves the multipole hierarchy and produces the theoretical predictions for the CMB temperature and polarization 
angular power spectrum, as well as the matter power spectrum. We have focused on the interaction Q£ = — rp c (l + 5 c )u^, 
where T is a co nstant rate of energy density transfer [see Eqs. (fl]) and Generalising a previous result for non-interacting 
dark energy in iDoran et al.l (|2003l ). we find that, in our interacting model, requiring adiabaticity between all the other 
constituents (photons, neutrinos, baryons, and cold dark matter) leads automatically also to dark energy adiabaticity, if its 



20081 ), 



early-time equation of state parameter is w e < —1 or —4/5 ^ w e ^ 1/3. In our previous work (jValiviita et al. 
showed that if the equation of state parameter for dark energy is —1 < Wde < —4/5 in the radiation or matter eras, the model 
suffers from a serious non-adiabatic instability on large scales. In this paper, the systematic derivation of initial conditions 
confirms that result. However, in this paper we have shown that the instability can easily be avoided, if we allow for suitably 
time- varying dark energy equ ation of state. The m ain results are verbally summarised in Table [1] o n page [5] 



mo 

In the companion paper ( Valiviita et al. 20091 ) we modified the CAMB Boltzmann integratoi^ ( Lewis et al. 2000h . 



using 

the adiabatic initial conditions derived here for the interacting model, and performed full Monte Carlo Markov Chain likelihood 
scans for this model as well as for the non-interac ting (r=0) model for a reference, with vari ous co mbinations of publicl y 
availa ble data sets: WMA P jKomatsu et al.ll2009l) . WMAP & ACBAR jReichardt et alj|2009h . SN jKowalski et al.ll2008l ), 
BAO jPercival et al]|2007l ), WMAP&SN, WMAP&BAO, WMAP&SN&BAO. 

With the parametrization Wde = woa + w e (l — a), viable interacting cosmologies result fo r wo close to —1 and w e < — 1 
or —4/5 < w e ^ 1/3, as long as Wo + 1 and w e + 1 have the same sign ( Valiviita et al. 20091 ). These particular conclusions 
apply exclusively to the interaction model we considered in this paper. 

However, the method can be easily adapted for studying different interactions: one only needs to modify the background 
evolution and interaction terms in Eqs. (|29p . (|30p . (|37p . and (|38p . before reading a new matrix A(a;) from them. Based on 



http:/ /camb.info 
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section IV of Valiviita et al. ( 20081 ) . the other interacting fluid models [aQ c = ~aTLp c or aQ c — —(3TL{p c + Pde), where 



a, /3 < 1 are dimensionless constants], that are common in the literature, behave in a very similar way to the model studied 
here, i.e., for —1 < w e < w C rit the models are not viable due to the early-time large-scale blow-up of perturbations, for 
Wcrit < We < wWiab the models can be viable and non-standard adiabatic initial conditions maybe found, and for w e > u> a diab 
(or We < — 1) the models are viable and standard (non-interacting) adiabatic initial conditions can be found. The critical value 
Wait is determined by demanding that the 'blow-up' mode is actually a decaying mode and the fastest 'growing' curvature 
perturbation mode i s a constant, i.e., that the largest real part of the eigenvalues Xi is zero, which with the notation of 



Valiviita et al. ( 2008h is guaranteed whenever Re(n+) ^ 3. In our model the critical value, iu C rit = —4/5, is independent of 



th e strength of in t eract ion, but in the above-mentioned models it depends on a or /3, as indicated by equations (85) and (98) 
in Valiviita et al. | |2008h . In general, our results show that the (early-time) dark energy equation of state plays, together with 



the interaction model, an important role in the (in)stability of perturbations. 
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